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Abstract. - We present molecular dynamics simulations on strongly charged polyelectrolytes 
in poor solvent. The resulting pearl-necklace conformations are analyzed in detail. Fluctuations 
in the number of pearls and their sizes lead only to small signatures in the form factor and the 
force-extension relation, which is a severe obstacle for experimental observations. We find that 
the position of the first peak in the structure factor varies with the monomer density as ~ Pm 35 
for all densities. This is a qualitative difference to polyelectrolyte solutions in good solvent 
which scale &s Pn{ and in the dilute and semi-dilute concentration regime, respectively. 



Polyelectrolytes (PEs) are polymers that carry ionizable groups that dissociate ions in 
aqueous solution. Technical applications range from viscosity modifiers, precipitation agents, 
superabsorbers to leak protectors [Jj] . In biochemistry and molecular biology they are of great 
importance because virtually all proteins, as well as DNA, are PEs. 

Many PEs contain a non-polar hydrocarbon backbone, for which water is a poor solvent. 
Therefore, in aqueous solution, there is a competition between the tendency to precipitate, the 
Coulomb interaction and the entropic degrees of freedom. This can lead to elongated strings 
of locally collapsed structures (pearls) . Such necklace conformations have been predicted on 
the basis of scaling arguments in ref. 0] for a weakly charged single chain PE, and have 
been confirmed by simulations using the Debye-Hiickel approximation and with explicit 
counterions jij. The formation of the necklace structure is due to the Rayleigh instability 
of a charged droplet, which leads to a split once a critical charge is reached. The size of 
the pearls is determined by the balance between electrostatic repulsion and surface tension. 
The distance between two pearls is governed by the balance of the electrostatic pearl-pearl 
repulsion and the surface tension. However, there is up to now no clear experimental proof for 
the existence of necklace chains Our previous papers dealt with smaller system sizes 

and showed that necklaces exist also in the presence of counterions and exhibit a variety of 
conformational transitions as a function of density. The focus of the present paper is to analyze 
by extensive computer simulations in detail three possible experimental observables, namely 
the form factor, the structure factor and the force-extension relation, which can be probed 
by scattering and AFM techniques. Analyzing the fluctuations of necklace structures we find 
an extended coexistence regime between different necklace structures and broad distributions 
for the pearl sizes and the pearl-pearl distances that smear out the necklace signatures. In 
addition we find that the peak in the structure factor of the solution scales proportional to 
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Fig. 1 - Snapshot of a necklace structure analyzed with the cluster algorithm. Chain length N m = 382: 
(black) monomers in pearls, (light gray) monomers in strings, (gray) counterions. 



1/3 

the density p n i from the dilute up into the dense phase which is in clear contrast to good 
solvent PEs. 

Our PE model is the same as in Ref. and consists of a bead spring chain of Lennard- 
Jones particles, whose interaction at distance r < R c is given by 4e[(-) 12 — (-) 6 — c], and is 
zero elsewhere. The constant c is chosen such that the potential value is zero at the cutoff 
i? c , and e is a measure of the solvent quality. Chain monomers interact up to a distance 
R c = 2.5(7 and experience an attraction with e = 1.75 k^T. The 0-point for this model is at 
e = 0.34fceT [Q. The counterions interact via a purely repulsive interaction with R c — 2 1 / 6 a\ 
The units of length, energy and time are a, e, and r, respectively. For bonded monomers 
we add a bond potential of the form — 14fcB?Tn(l — (^) 2 ); which results in an average bond 
length b = 1.09a. Charged particles at separation ry = |rj — v\\ interact via the Coulomb 
energy fceT '^B<Mj / 't %j , with qi = 1, (—1) for the charged chain monomers (counterions). The 
Bjcrrum length £b = e 2 / (AweseokBT) (e: unit charge, eo and es: permittivity of the vacuum 
and of the solvent) is set to £b = 1.5er. A velocity Verlet algorithm with a standard Langevin 
thermostat is used to integrate the equation of motion @ (friction coefficient T — t~ x , time 
step At — 0.0125t). The Coulomb interaction was calculated with the P3M-algorithm ||, 
tuned to force accuracies well above the thermal noise level. The central simulation box of 
length L contained Np — 5, 10 or 32 chains subject to periodic boundary conditions. Each 
chain consisted of N m = 48 . . . 478 monomers, with a charge fraction / = 1/3 ... 1/2. Other 
parameter sets yield qualitatively similar results ^(Tof, so that we restrict ourselves here to 
the above described "generic" system. The simulation time for the longest chains was around 
200 times the measured correlation time for the end-to-end distance Re and the chains centers 
of mass diffused at least several radii of gyration Rq- The pressure p is always positive and 
the pV diagram is convex at all densities, thus our simulations are stable, reach true thermal 
equilibrium, and reside in a one phase region. The volume density inside the pearls does not 
exceed 0.47, which is well below the glass transition. 

Being first interested in single chain properties we used a monomer density p m — Np I ^ m = 
1.5 x 10~ 5 (T~ 3 . A typically 4-pearl (and 3-string) structure is shown in fig. |l|. The substructures 
contain 90-8-94-6-77-9-98 monomers. For the identification of the pearls and strings we used 
an especially adapted cluster algorithm For chains of length N m = 430, we find 

structures with 4, 5, 6 and 7 pearls and an average pearl number np = 5.22. Results for other 
chain lengths can be found in table |[ The coexistence regime between different structure types 
enlarges with increasing chain length. One reason for the small differences in the free energy is 
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Table I - Basic observables for various chain lengths at p m — 1.5 x 10~ 5 a -3 . Symbols: Chain length 
N m , end-to-end distance Re, average number of pearls np, distance between neighboring pearls rpp, 
number of monomers gp in a pearl, pressure p and the osmotic coefficient II. Types denotes the 
observed structures with np pearls, and Srpp and Sgp denote the relative standard deviations. All 
systems had a Bjerrum-length £g = 1.5a and a charge fraction f = 1/3. 
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the interplay between the chain conformation and the counterion distribution. Conformations 
with a lower number of pearls have a smaller extension, e.g. R B , and a larger pearl size. 
This leads to a stronger attraction of counterions and yields a smaller effective charge on the 
chain which in turn stabilizes larger pearl sizes and smaller chain extensions. In contrast 
to scaling theories [pT[]l4]] we do not find a collapse into a globular state due to this effect. 
This is consistent with a more refined theoretical analysis jl5| that takes prefactors and finite 
concentrations into account. The interplay between counterion distribution and structural 
changes of the conformation can be measured e.g. via the number of counterions n c (n) for 
a n pearl structure inside a shell of 3a from the chain. For example, for N m = 430 we 
find n e (4) = 60.4, n c (5) = 56.7, n c {6) = 52.8 and n c {7) = 50.2. From the probability p(n) of 
finding a n pearl structure we calculate the free energy difference AT nm — kp,T In (p(n) / p(m)) . 
For N m = 430 we find AT 45 = -1.33/cbT, AT 56 = 0Mk B T, AT 67 = 1.90fc B T. All values 
are of the order k B T which is consistent with the observed large coexistence regime. In order 
to exclude the possibility of metastable states, we have observed the time evolution for a 
single chain. We find typically many transitions between different structure types, see i.e. one 
example for N m = 430 in fig. Not only the pearl number, also the position and size of the 
pearls fluctuate strongly. For chains with more than 2 pearls, both the mean values as well 
as the relative standard deviations 6 of the distance rpp between neighboring pearls and the 
number gp of monomers in a pearl are almost constant, see table |, which is consistent with 
scaling theory. 

Scattering experiments probe the conformation via the chain form factor Si (q) given by 

where (...) denotes the ensemble average. Fig. |] shows the measured form factor for 7V m = 382 
(see table |). In the Guinnier regime (Rcq "C 1) the radius of gyration Rq can be calculated 
from Si(q) = N m (l — (Rcq) 2 /3), giving Rq = 16.8<r ± 0.3cr, in agreement with the directly 
calculated value Rq — 16.9a ± 0.4er. In the range 0.07cr _1 < q < 0.3cr _1 Si scales as 
g~ - 93 which corresponds to a stretched object. Around q — 0.46cr _1 the form factor has 
an inflection point. A comparison with the intra pearl scattering reveals that this is due to 
inter pearl scattering (inset in fig. |^). Dividing out the intra pearl form factor gives access 
to the inter pearl scattering and thus rpp. The form factor yields rpp = 13.6er, in accord 
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Fig. 2 - Time evolution of a single chain of length N m = 430 in equilibrium, shown over a quarter of 
the physical running time. Many transitions between different structure types can be observed. 

with the directly measured value rpp = 13.3cr (table |). In the region around q = l.Ocr -1 we 
find Si(q) oc q~ 4 , the typical Porod scattering. From the small dip at q = \.7o~ x one can 
calculate the radius rp of the pearls to be rp = 2.6c which again compares well to the real 
space value rp ~ 3a. We conclude that the cooperative effect of fluctuations on overlapping 
length scales broadens all characteristic signatures which can be revealed by scattering under 
experimental conditions (polydispersity, charge fluctuations, etc). Thus necklaces might be 
difficult to detect. 

The overall scattering function S(q) of the solution contains additional experimental in- 
formation. We analyze here the inter chain scattering Sic — S/Si. For good solvent PEs 
experiments Jl6| , theory and simulations [i~8| find a pronounced first peak of Sic a t 

1/3 

q* = (27r)/£, where £ is the correlation length. The position varies as q* oc pm in the very 




Fig. 3 - Form-factor Si (q) for pearl necklaces: (solid) form-factor of the whole chain, (dashed) form- 
factor due to intra pearl scattering, (dotted) fits. The marked region is enlarged in the inset. 
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Fig. 4 - Density dependence of the peak q* in the structure factor for three different chain length 
N m = 100, 200, 300 with / = 0.5. The black line is a fit to the data with N m = 200. 



dilute regime and crosses over to a p]l 2 regime at higher concentrations. In fig. [| we have plot- 
ted the density dependence of q* in poor solvent for different chain lengths. Within the error 
bars we find that for poor solvent chains q* scales proportional to p m 35±0 - 04 for all concentra- 
tions and chain lengths. Possible reasons for this are that the response of the polyelectrolyte 
conformation to density changes is much larger in the poor solvent case than in the good 
solvent case |Q , and the chain extension behaves non-monotonic as a function of density . 
Furthermore, in the density regime between p m = 1CP 2 <7 -3 . . . 10~ 4 er~ 3 the chain extension 
and the pearl number varies most, and almost all monomers are located within the pearls. 
Upon approaching the dense regime, the string length tends to zero and we find a chain of 
touching pearls, indicating that the conventional necklace picture breaks down. Our result is 
compatible to scaling exponents found in scattering experiments p9[ . 

Scaling theories jl2],|2(J] have predicted a p][ 2 regime to start at p*, which is defined at 

1 /3 

the density where Re ~ £, and to extend until £ « rpp where a bead-controlled p n { regime 
starts. We find p* Q ~ 5 x 1CT 2 , 1CT 3 , 1CT 4 for N m = 100,200,300, and r PP is equal to £ 
between p m = 10~ 2 and 10 _1 . One probable reason for our different findings is the strong 
inter-chain coupling and influence of the counterions on the conformations which are not 
sufficiently taken into account in the mean field approach. It is not clear at this stage if the 

1/2 

Pm regime can be recovered for much longer chain length. In addition we observe that the 
chains form a transient physical network at p m = 0.2er~ 3 for iV m > 200 which has neither 
been seen in previous simulations nor predicted by theoretical approaches but is in accord 
with experimental studies ||l9| . During the simulation time these networks reconstruct several 
times, e.g. chains are not trapped. 

A different way to measure a necklace signature can be obtained by single-molecule force 
spectroscopy, e.g. stretching the chain by AFM plj or optical tweezers. The case of an 
imposed end-to-end distance was investigated in Refs. JT3|,p2],|2^] for a weakly charged chain 
at infinite dilution, for which a saw-tooth pattern in the force-extension curve was predicted. 
We performed simulations in the weak coupling limit, using a single chain with N m = 256, 
/ = 1, £b = 0.08ct, and no counterions present. At equilibrium this system is in a two-pearl 
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Fig. 5 - Force-extension relation for a transition from two to three pearls. Shown are the average 
over all conformations, over only two-pearl configurations and only three-pearl configurations. 



configuration with Re = 21a. The force-extension curve was obtained by imposing various 
fixed i?E up to Re = 50a. Then the force on the end-monomers was measured, see fig. |^. The 
first remarkable observation is that the two-pearl state evolves into a three-pearl state under 
extension, which was predicted in ref. [^2[|2^]. This is counter intuitive from simple arguments, 
but can be shown to result from the electrostatic inter-pearl interactions Only when the 
force averages where computed separately for the three- and two-pearl states, we recognize the 
predicted saw-tooth pattern in the force. In equilibrium one would expect a rounded plateau 
for the transition, and the upper and lower part would correspond to metastable superheated 
and supercooled states. Again, due to the large fluctuations all nontrivial signatures in the 
force-extension relation are washed out. Analyzing the variance of the measured forces by 
applying a time window of the order of the transition time we find a 25% increase in the 
width of the force distribution in the transition region Re = 30 . . . 34c, and a bimodal force 
distribution. This means that one needs an experimental time resolution below one us to 
directly resolve the structure dependent force difference. An estimate of typical viscous and 
inertial forces acting on colloidal particles shows that the Brownian motion of small spheres 
attached to a single PE chain might indeed reflect the fluctuating forces exerted by the chains. 
The experimental determination of theses fluctuating forces by optical nanorheology 0] or 
AFM noise analysis |25|] presumably is difficult due to the high characteristic frequencies 
involved. 

To sum up we find that poor solvent PEs show an extended coexistence regime between 
different necklace structures. We observe stable necklaces even with condensed counterions. In 
all investigated cases the signatures of necklaces in the chain form factor are broadly smeared 
out. A similar broadening of signatures is found for the force-extension curve which will 
obscure any of the recently predicted saw-tooth patterns [|l3],^2],^j|. Nevertheless, the simu- 
lations reproduce the predicted initial increase of pearl number under imposing an external 
length constraint |^2|,|2^]. Our results suggest that the strong fluctuations impose a severe 
obstacle in observing necklace structures and structural transitions directly in experiments. 
Finally the location of the first peak in the solution structure factor is qualitatively different 
from that of PEs in good solvent. This is in good accordance to experiments jli| but does 
not show the full complexity of the proposed scaling picture at this point [ po[ . 
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Note: After the submission of this work we learnt that recent light scattering measurements 
of D. Baigl and C. E. Williams of q*(p m ) for PSS solutions produce an exponent of 0.38 for 
a charge fraction of 55% which is in good agreement with our simulations. 
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